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In strongly correlated electron systems, the emergence of states in the Mott gap in the single¬ 
particle spectrum following the doping of the Mott insulator is a remarkable feature that cannot 
be explained in a conventional rigid-band picture. Here, based on an analysis of the quantum 
numbers and the overlaps of relevant states, as well as through a demonstration using the ladder 
and bilayer t-J models, it is shown that in a continuous Mott transition due to hole doping, the 
magnetically excited states of the Mott insulator generally emerge in the electron-addition spectrum 
with the dispersion relation shifted by the Fermi momentum in the momentum region where the lower 
Hubbard band is not completely filled. This implies that the dispersion relation of a free-electron-like 
mode in the electron-addition spectrum eventually transforms into essentially the momentum-shifted 
magnetic dispersion relation of the Mott insulator, while its spectral weight gradually disappears 
toward the Mott transition. This feature reflects the spin-charge separation of the Mott insulator. 

PACS numbers: 71.30.-l-h, 71.10.Fd, 74.72.Gh, 79.60.-i 


I. INTRODUCTION 

A significant feature in the Mott transition is the 
emergence of states in the Mott gap in the single-particle 
spectrum due to the doping of a Mott insulator/i Since 
such a feature does not appear when a noninteracting 
band insulator is doped, this phenomenon must be 
related to strong electronic correlations intrinsic in 
the Mott transition. However, although these states, 
which are called doping-induced states or in-gap states, 
have been observed in materials such as cuprate high- 
temperature (high-Tc) superconductors and in various 
theoretical calculationstheir interpretations are 
controversial. 

In particular, the question of whether the doping- 
induced states in the two-dimensional (2D) Hubbard 
and t-J models in the parameter regime relevant to 
high-Tc cuprates near the Mott transition are essentially 
disconnected from the low-energy states by an energy 
gap or not has attracted considerable interest. This is 
because this issue is related not only to the nature of 
the Mott transition but also to the Fermi surface near 
the Mott transition (hole pockets or a free-electron-like 
Fermi surface)^d^ ~^^'^"^~^^'^^i^^i^^~ — and the anomalous 
features observed in high-Tc cupratesA"— The presence 
of the energy gap has been suggested by numerical 
calculations^i"— i22ii22r— and interpreted in terms of 
various concepts, such as the binding between double 
occupancy and vacancy^i^i^ii^iSl the effect of a charge 
2e bosonj^ii^ hybridization between the quasiparticle 
and cofermion)^i21 and a spin-polaron shakeoffi^i^ The 
absence of the energy gap has also been suggested by 
numerical calculations ) where the doping- 
induced states have been interpreted as a quasiparticle 
coupled with antiferromagnetic correlationsi^ and as 
essentially the magnetic excitation of the Mott insulator 
which emerges in the single-particle spectrum because 
charge character is added by dopingj^^— 

In this paper, we address this issue from a general 


perspective, by investigating how the dispersion relation 
of the doping-induced states should generally behave in 
the small-doping limit of a continuous Mott transition, 
based on an analysis of the quantum numbers and the 
overlaps of relevant states. The validity of the argument 
is demonstrated using the two-leg ladder and bilayer t-J 
models. An intuitive real-space picture of the doping- 
induced states is also presented. 

We consider zero-temperature properties of repulsively 
interacting systems where the ground states are Mott 
insulators at half-filling and metals in the small-doping 
limit. The numbers of sites, unit cells, electrons, and 
doped holes are denoted by Ng, N^, and iVh(= Ng — 
Ne), respectively. The doping concentration is denoted 
by 6{= Ni,/Ng). 

II. QUANTUM NUMBERS 

We study the single-particle spectral function defined 
as 

A(k oj]=l 5Em,^K’^l4,alGS)/*P<5(a;-e™) for w > 0, 

’ \ lJ2m,a\iM(^kMGS)h\‘^S{uJ 4-Em) for W < 0, 

where |GS)ft, and Sm denote the ground state at 
fVh = h and the excitation energy of the eigenstate 
|m) from |GS)?i, respectively. Here, c], ^ denotes the 
creation operator of an electron with momentum k and 
magnetization a. 

If |GS)?t has spin s, magnetization Sz, and momentum 
fcp, ^|GS)/i can be nonzero only for \m) with 

h — 1 holes, spin |s ± 1/21, magnetization Sz + cr, and 
momentum fc -f In addition, if the ground-state 

energy (including the chemical-potential term) at Ni, = h 
is lower than that of A^h = h — 1 by e, iVe-conserving 
excited states with excitation energy ui — e from |GS)/i_i 
can contribute to A{k,uj) at Ah = h, where e —>■ -1-0 
for As —>■ oo in a metallic phase (including h = 1). 
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In interacting systems, it is generally expected that 
low-energy eigenstates (particularly those of dominant 
excitations) having the same quantum numbers as those 
of components of cj[^|GS)ft, have nonzero overlap with 

^|GS)/j for k where the lower Hubbard band (LHB) is 
not completely filled. 

Thus, if |GS)i has spin 1/2 and momentum fcp, |GS)o 
with spin 0 and momentum 0 should contribute to 
A{k,uj) at fe = —fcp and w = e for the top of the hole¬ 
like band [LHB at half-filling (A^h = 0)] in the small¬ 
doping limit, which causes hole-pocket-like behavior as in 
a doped band insulator)^ In addition, the magnetically 
excited states with spin 1 at half-filling^ having the 
dispersion relation w = /(fe) should emerge following 
doping along w = /(fe -I- fep) -b e in A{k,oj) for fe where 
the LHB is not completely filled. 

According to this argument, if the magnetic excitation 
of the Mott insulator is gapless at momenta 0 and Q, the 
doping-induced states in the small-doping limit should 
also be gapless at —fep and Q — kp (±fep if Q = 2fep) 
and exhibit the momentum-shifted magnetic dispersion 
relation in the momentum region where the LHB is not 
completely filled (primarily outside the Fermi surface). 
This behavior essentially agrees with the results for the 
one-dimensional (ID) and 2D Hubbard models in Refs. 

and those for the 2D t-J model in Ref. 
and suggests the coexistence of hole-like behavior and 
electron-like behavior for w > 0. 

Although it has been pointed out that the spin-0 and 
spin-1 states at half-filling can appear in the electron- 
addition spectrum at Np = 1, the relationship of the 
momenta (dispersion relations) has not been recognized 
in Ref. [13; the doping-induced states in the 2D t-J and 
Hubbard models have been considered to be disconnected 
by an energy gap from the low-energy states relevant 
to hole-like Fermi surfaces and have been interpreted as 
spin-polaron shakeoff bands in Refs. iaii However, 
according to the above argument, it is natural to interpret 
the doping-induced states in the small-doping limit as 
essentially the states of the spin-wave mode of the Mott 
insulator and to consider that they exhibit the gapless 
spin-wave dispersion relation^^— shifted by fep in the 
2D t-J and Hubbard models. 


III. OVERLAP 

To strengthen the above argument, we consider the t-J 
model defined by the following Hamiltonian: 

Ti = 'y ( + y ^ Ji,j{Si ■ Sj — —niTlj) — fj, ^ ( Tli, 

ijj i 

where = Ci,o.(l - n^-a) and nt = for the 

annihilation operator Ci^a and number operator of 
an electron with magnetization a at site i. Here, Si 
denotes the spin operator at site i. At each site, a vacant 
state jo) or a singly-occupied magnetization-cr state jcr) 


is allowed, but double occupancy is forbidden. Here, we 
assume that |GS)o has spin 0 and momentum 0 and that 
|GS)i has spin 1/2, magnetization and momentum 
fep. For instance, |GS)o on a bipartite lattice with the 
same number of sites in the two sublattices has spin 0 
in a finite-size system with Jij > 0 for neighboring sites 
between the two sublattices and Jij = 0 otherwisei^ 
Because cj^|a)* can be expressed as _^\aY and 

using 


r 


qZ,(T 

,-CT 5 




where jo:)* denotes the state at site i (c|ctI“)* = 1^)* 
\aY = |0)/ and 0 otherwise). 


o(GS|cf_fe^__JGS)i 


where 


:^Sjo(GS|5' 


— kp—ki 

— kp—ki 


4., JGS)i 

4., -JGS)i, 

( 1 ) 


‘"k.a 








using the site j coordinate rj. Equation (1) means that 
the overlaps between the electron-addition states at Np = 
1 (cj^^|GS)i and _^|GS)i) and the single-spin excited 

states at half-filling (S'^/fcp|GS)o and |GS)o) are 

related to the overlap between the electron-removal state 
for the momentum at the top of the LHB at half-filling 
(c_fc,,_,|GS)o) and |GS)i. 

From Eq. ([T]), if o(GS|c[Lfcp__^|GS)i is 0(1), 
o(GS|5X-,4_JGS)i and o(GS|5!’,7_,4__JGS)i are 
expected to be 0(1/'/Ns) at 0{Ns) fe-points. Then, for 
these fe-points, the normalized electron-addition states 
should normally have 0(1) overlap with the normalized 
single-spin excited states because o(GS|5'„fc. • 

Sk,\GS)o = 3/4 and E, i(GS|dfc...4.,JGS)i = 
1/iVs^d^ 

In addition, because all sites are singly occupied 
at half-filling, cJ,< 7 G>' |GS)o = |GS)o and 

= 4,fcp|GS)o. Hence, 


oiGSlSX-kcljGSY 

o(GS|.5E7-.4.-JGS)i 


VKo{GS\cX^_^CkAjG^)i, 

VKo{gs\cX^_X.X,-Xsy 

Sk,-k^VKo{GS\cX -JGS)i/2. 


Thus, the overlap can be roughly estimated in a 
decoupling approximation as 

o{GS\sYC-A,YGSY « v^ei(GS|c4,.4jGS)i, 
o(GS|^E“-fc4.-JGS)i « v^ei(GS|c4,-.4_JGS)i 

( 2 ) 

where ^ = o(GS|c[Lfcp _^|GS)i. At fe = —fep. 


o{GS\sYX-AJGS), = o(GS|.SE”-fe4,-JGS)i = 0 
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because |GS)o is assumed to be a spin-0 state. 

The above results imply that the single-spin excited 
states of the Mott insulator contribute significantly to 
the doping-induced states in the small-doping limit for 
k where the LHB is not completely filled (cj[ ^ | GS) i ^ 
0) and for k ^ -fcp if o(GS|cL^,j, _JGS)i is 0(1). 
In ID systems, even if o(GS|c^^p _^|GS)i —>■ 0 for 
TVs —>■ 00 ^ 1 ^ exact-solution analyses and numerical 
calculations strongly suggest that the mode of the 
doping-induced states in the small-doping limit primarily 
reflects the dominant part of the magnetic excitation.— 


IV. WEAKLY COUPLED TWO-SITE 
CLUSTERS 

To demonstrate the validity of this argument, we 
consider the two-leg ladder and square-lattice bilayer t-J 
models for ty, Jj_ Jy, tj. J_l, and ty Jy > 0 
(J_L/ty is not too small for the ground state to have 
spin 0 or 1/2 [Hii). Here, t± and Jl denote —tij 
and 2Jij for neighboring sites between chains (planes), 
respectively, and ty and Jy denote —Uj and 2Jij for 
neighboring sites in chains (planes), respectively, in the 
ladder (bilayer) model. The other tij and Jij are set to 
zero. Although some aspects of the doping-induced states 
in similar models have been discussedj^^— this paper 
clarifies the nature of the doping-induced states through 
explicit calculations of the spectral weights, dispersion 
relation, and overlaps with magnetically excited states, 
using the effective eigenstates up to 0(ty) and 0(Jy). 


A. One-hole doping 


The ground state at half-filling (spin 0 and momentum 
0 ), the low-energy spin-1 eigenstate at half-filling, and 
the low-energy eigenstate at iVh = 1 are expressed, 
respectively, as 

|GS)o = |S)"“, 


where [S)-^ , |Ts^)l, and represent the two-electron 

spin-0 state, spin-1 state, and one-electron bonding 
state, respectively, at the j-th two-site unit cell with 
coordinate ryj for ty = Jy = 0. Here, Sz and fey 
denote the magnetization and inter-unit-cell momentum, 
respectively. The magnetic dispersion relation at half¬ 
filling is obtained by calculating the excitation energy of 
|T)o"’*^ from |GS)o: 


AA(fey) = | 


Jy cos fey -j- J_i_ 

Jy (cos feya; ~\~ cos fey^) Jj_ 


for ID systems, 
for 2D systems. 


The ground state at Vh = 1 with magnetization (|GS)i) 
is if fey can be set to tt, where tt indicates tt 


and (tt, tt) for ID and 2D systems, respectively. For 
the small-doping limit, the value of /r is set such that 
i(GS|'H|GS)i = o(GS|'H|GS)o - e with e ^ -HO. 

We calculate the spectral function at iVh = 1. By 
adding an electron to the ground state at Yh = 1, spin- 
0 and 1 states are obtained^ as follows: |GS)o/\/2(= 

cLfcp_JGS)i) for spin 0, and |T)p“^’'’^Vv^(= 
£(fe„..),JGS)i) and |T)o""+"’VV2fV;:(= 
for spin 1. Here, (fey,fe_L) indicates the momentum 
with inter- and intra-unit-cell momenta fey and fej_, 
respectively, and fep = (7r,0). Thus, there are 

spectral weights at fey = —tt and w = e in Ab(fey,w) 

lo(GS|cj[ ,^|GS)ip = 4-fcp/4] and essentially 
along the dispersion relation of the magnetic excitation 
shifted by TT, a; = A£'(fey -H tt) -|- e, in Aa(fey,w) 

[iEs.,^lo"’^’'’""(T|S(fc||,^),<,|GS)ip = 3/(4Vu)], where 
Ab(fey, w) and Aa(fey, w) denote the spectral functions for 
^(fc|i 0 ) <T tt) ct’ respectively. The shift of tt is due 

to the momentum difference between |GS)i and |GS)o. 

The spectral weight for w > 0 averaged over the 
momentum is 1/W(= ^), which is consistent with the 
spectral-weight sum rule— ^ One quarter of this spectral 
weight (5/4) is from the ground state at half-filling, 
which contributes to the mode causing hole-pocket-like 
behavior, whereas three quarters (35/4) is from the 

magnetically excited states (|T)q"’^^), which spread over 
the entire fey region for fej_ = tt and can be identified as 
the doping-induced states. In addition, Eq. o and the 
validity of the decoupling approximation [Eq. ([2])] can be 
confirmed using 


o(GS|cLfc, _JGS)i 


l/A 

6k^,^/V^s. 


B. Multi-hole doping 

A similar analysis is applicable to multi-hole-doped 
cases in the small-doping regime. The ground state at 
Ah = fe is effectively expressed as a linear combination of 
direct products of h jB^) and {N^ — h) |S) states. Because 
the doping-induced states are obtained by replacing one 
of the h |Bcr) states with jT^^), the spectral weights 
increase as 5 increases. In addition, because a magnetic- 
excitation mode similar to that of the Mott insulator 
remains for a; ~ J_l and because the Fermi momenta 
remain near (7r,0), significant spectral weights should 
continue to be located almost along the momentum- 
shifted magnetic dispersion relation of the Mott insulator 
in Aa(fey, to). The states obtained by replacing one of the 
h |Bct) states with |S) form a mode as in a doped band 
insulator (Ri the upper edge of the LHB at half-filling) in 
Ab(fey,a;)4a 

More simply and generally, the behavior of the spectral 
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FIG. 1: (a) Aa(fc||,‘^)^|| at A^h = 2, (b) ^a(fe||,at A^h = 0, and (c) 5a(fc||,a;)t|| at A^h = 0 in the t-J ladder for tx/t|| = 2, 
J||/t|l = 1/4, Jx/t|| = 1, and Ns = 120, obtained using the non-Abelian DDMRG method. At A^h = 0, /r is set to the value of 
the small-doping limit, (d)-(f) The same as (a)-(c) but for Ab(fc||,a;)t|| and Sb(fc||,(^)t||. The solid green lines indicate uu = 0. 
Gaussian broadening with a standard deviation of 0.1f|| is used. 


function at finite S{= So) can be understood as being 
deformed continuously from that of the small-doping 
limit if phase transition does not occur for 0 < (5 < ^o- 
This change should be small in the small-doping regime. 

For confirmation. Fig. [T] shows the results for the 
t-J ladder obtained using the non-Abelian dynamical 
density-matrix renormalization group (DDMRG) 
method^ [the DDMRG method^ in the U(1)G>SU(2) 
basi o^’^'^^ ], where 120 density-matrix eigenstates are 
retained. The spin dynamical structure factors S'a(fc||,a;) 
and 5'b(fc||,a;) are defined, respectively, in the same 
way as Aa(fc||,a;) and Ab(fc||,w) for w > 0 but for spin 
operators. The mode for w/t|| « 1 in Aa(fc||,a;) at 
A^h = 2 corresponds to the doping-induced states and 
exhibits almost the same dispersion relation as that 
of the magnetic excitation at A^h = 0, but shifted by 
approximately tt in [Figs. [Il[a)-[T][c)]. The states for 
fc|| Ri TT and w Ri 0 in Ab(fc||,a;) form a mode as in a 
doped band insulator [Figs. [T][d) and[T][e)]. Thus, the 
results shown in Fig. [T] are consistent with the above 
analysis. 

The results shown in this section fSec. lIVI) demonstrate 
that the magnetically excited states of the Mott insulator 
emerge as the doping-induced states in the small-doping 
limit in ID and 2D systems and that the dispersion 
relation of the doping-induced states is essentially the 
same as that of the magnetic excitation, but shifted by 
the Fermi momentum. In the above ladder and bilayer t- 
J models, the doping-induced states in the small-doping 
limit are disconnected from the low-energy states because 
the magnetic excitation at half-filling has an energy gap; 
in models having a gapless magnetic excitation at half¬ 
filling, the doping-induced states should exhibit a gapless 


dispersion relation in the small-doping limit. 

V. DIFFERENCE DUE TO DOPING 

The difference in the electron-addition spectra between 
a Mott insulator and the small-doping limit is explained 
as follows. In a Mott insulator, c^^|GS)o yields no 
spectral weights in the Mott gap because there are no 
eigenstates with A^e = -^s + 1 in this energy regime. 
In the small-doping limit, in contrast, cj[^|GS)i can 
yield spectral weights in the Mott gap because low- 
energy eigenstates for the magnetic excitation of the 
Mott insulator exist 

The emergence of the electron-addition states in the 
Mott gap reflects the spin-charge separation of the 
Mott insulatori^^— The spin-charge separation means 
that the low-energy properties are described in terms 
of independent spin and/or charge degrees of freedom 
rather than in terms of an electron-like quasiparticle. 
In ID interacting electron systems, the spin-charge 
separation occurs even in a metallic phase in the low- 
energy limit<^“— On the other hand, in a Mott insulator, 
the spin-charge separation occurs^ in any dimension: 
the low-energy properties are described in terms of the 
spin excitation which is gapless or has an energy gap of 
the order of the (effective) spin-spin interaction, whereas 
the charge excitation has an energy gap of the order of 
the Coulomb interaction. Because of this spin-charge 
separation, although there are low-energy magnetically 
excited states, the electron-addition excitation does not 
have spectral weights in the Mott gap in a Mott insulator; 
the charge quantum number of the magnetically excited 
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states (iVe = Ns) differs from that of cj[^|GS)o {N^ = 
Ns + 1). However, the magnetically excited states 
can emerge in the electron-addition spectrum with the 
dispersion relation shifted by the Fermi momentum 
following the doping of the Mott insulator as discussed 
abovei^^— This is in contrast with the noninteracting 
band insulator case, where the lowest spin excitation 
energy is the same as the band gap. 

In the ladder and bilayer t-J models considered above, 
although the ground-state behavior is similar to that 
of a (doped) band insulator, the characteristic of the 
Mott transition (the emergence of the states) occurs for 
a; > 0, reflecting the spin-charge separation of the Mott 
insulator (the presence of the low-energy spin excitation 
and the absence of the low-energy charge excitation). 

VI. REAL-SPACE PICTURE 

The situation in the small-doping regime can be 
intuitively explained as follows. The added electron, 
which is located at the position of a vacancy at low 
energies, has almost no chance to hop to another vacant 
site if vacancies are distributed homogeneously. However, 
the added electron can hop with a small energy to 
a neighboring site occupied by an electron having the 
opposite magnetization, at the same time as that electron 
hops to the site that was occupied by the added electron. 
This is nothing but the spin exchange process. Thus, 
in general, the electron-addition excitation in the small¬ 
doping regime exhibits a dispersion relation similar to 
that of the magnetic excitation of the Mott insulator, but 
shifted by the Fermi momentum. This shift is due to the 
momentum difference between the ground states before 
and after the electron addition. The spectral weight in 
the large-repulsion regime should be 0{6) because the 
probability that the added electron occupies a vacant site 
is essentially Ni^/Ns in the small-doping regime. 

VII. REMARKS 

Even if the ground state at half-filling has an 
antiferromagnetic long-range order, c][£,.|GS)i should 
have nonzero overlap with the spin-wave mode at 
half-filling because both involve spin-1 components if 
|GS)o and |GS)i have spin-0 and spin-1/2 components, 
respectively. Thus, the doping-induced states in the 
small-doping limit are expected to exhibit essentially the 
same dispersion relation as that of the spin-wave mode, 
but shifted by the Fermi momentum, in the momentum 
region where the LHB is not completely filled (primarily 
outside the Fermi surface). This picture is consistent 
with the results for the 2D Hubbard and t-J models in 
Refs. IMM 

The continuous transformation of the single-particle 
dispersion relation for w > 0 into the momentum-shifted 
magnetic dispersion relation of the Mott insulator does 


not necessarily mean that the magnetic excitation of 
the Mott insulator is interpreted in terms of spinons 
(spin-1/2 and charge-0 quasiparticles), because this 
characteristic is independent of quasiparticle pictures. 
The momentum shift is simply due to the momentum 
difference between the ground states before and after the 
electron addition. In the above ladder and bilayer t-J 
models, the magnetic excitation is typically interpreted 
in terms of a triplon (spin-1 and charge-0 quasiparticle). 

The doping-induced states in the small-doping limit 
can involve not only the magnetically excited states 
but also other low-energy excited states (spin-0 excited 
states) at half-filling if cj[,^|GS)i contains components 
with these quantum numbers. Nevertheless, the doping- 
induced states are expected to be primarily due to the 
magnetic excitation because this is generally the most 
dominant low-energy excitation in a Mott insulator (Fig. 

In a first-order Mott transition where phase separation 
occurs for 0 < (5 < 6c{j^ 0), the states obtained by 
adding an electron in the 5 —>■ 0 limit are mixtures of 
the electron-addition states from the ground state at 
6 = Sc and those at i5 = 0, which are not exactly the 
magnetically excited states at half-filling. Nevertheless, 
the dispersing mode for w > 0 in the LHB before 
the Mott transition in the large-repulsion regime should 
more or less reflect the magnetic excitation from a Mott 
insulating state (an eigenstate at half-filling) similar to 
the ground state at (5 = Sc- Thus, the indications 
that a spin-charge-separated Mott insulator is being 
approached, such as a reduction in spectral weight 
from the dispersing mode,— are expected to generally 
appear near the Mott transition in the large-repulsion 
regime, provided the value of Sc is relatively small. 


VIII. DISCUSSION AND SUMMARY 

We consider the issue of the doping-induced states 
in the 2D Hubbard and t-J models mentioned in 
Sec. D Because the low-energy magnetic excitation 
at half-filling is well described in terms of the gapless 
spin-wave mode^^^— according to the present results, 
a mode exhibiting the spin-wave dispersion relation 
shifted by the Fermi momenta should emerge in the 
electron-addition spectrum in the small-doping limit; the 
mode should be gapless at the Fermi momenta because 
the spin-wave mode is gapless. If the LHB around 
(0,0) is essentially completely filled, the mode appears 
primarily outside this momentum region. This picture is 
consistent with the results in Refs. I33l - l35l : interpretations 
suggesting that the doping-induced states are essentially 
disconnected from the low-energy states by an energy 
gap2i“— i2Zr— are not consistent with this picture in the 
small-doping limit. For the ID Hubbard model, this 
picture essentially agrees with the results in Refs. 1^1^ . 

In summary, the general and direct relationships of 
the doping-induced states with the magnetically excited 
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states of the Mott insulator were clarified based on 
an analysis of the quantum numbers and the overlaps, 
as well as through a demonstration using the ladder 
and bilayer t-J models. The results imply that the 
doping-induced states in the small-doping limit exhibit 
essentially the same dispersion relation as that of the 
magnetic excitation of the Mott insulator, but shifted 
by the Fermi momentum, in the momentum region 
where the LHB is not completely filled in a continuous 
Mott transition. This characteristic should be general 
and independent of the dimensionality and quasiparticle 
pictures. 

With the reverse use of this characteristic, the 
magnetic dispersion relation of the Mott insulator can 
be inferred from the dispersion relation of the electron- 
addition excitation in the small-doping limit by shifting 


this dispersion relation by the Fermi momentum. In 
particular, if the electron-addition excitation has a gap 
in an energy regime, the magnetic excitation is also 
expected to have a gap in this energy regime. 
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